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We reconsider the influence of two-particle final state in- 
teractions (FSI) on two-particle Bose-Einstein interferometry. 
We concentrate in particular on the problem of particle emis- 
sion at different times. Assuming chaoticity of the source, 
we derive a new general expression for the symmetrized two- 
particle cross section. We discuss the approximations needed 
to derive from the general result the Koonin-Pratt formula. 
Introducing a less stringent version of the so-called smooth- 
ness approximation we also derive a more accurate formula. It 
can be implemented into classical event generators and allows 
to calculate FSI corrected two-particle correlation functions 
via modified Bose-Einstein "weights". 



I. INTRODUCTION 

Two-particle correlations in momentum space can be 
used to extract information about the space-time struc- 
ture of the emitting source 0J|]. The method exploits 
in an essential way the quantum mechanical uncertainty 
relation between coordinates and momenta, and thus 
any formal treatment of two-particle correlations must 
be based on a quantum mechanical description. For so- 
called " chaotic" sources where the two particles are emit- 
ted independently the description can be based on the 
single-particle Wigner density S(x, K) of the source. It 
is, however, known to be important (at least in principle, 
but for sufficiently small sources also in practice) that 
one starts from the correct quantum mechanical Wigner 
density rather than directly from a classical space distri- 
bution S c -[ 3SS {x, K) of the particles in the source because 
the latter can lead to unphysical behavior of the correla- 
tion function jj). The correlations are "generated" after 
the emission by two classes of effects: 

(1) For pairs of identical bosons (fermions) the two- 
particle wave function describing their propagation to- 
wards the detector must be symmetrized (antisymme- 
trized). For boson pairs this results in quantum statisti- 
cal "Bose-Einstein correlations" between the final state 
momenta of the two particles. Via the uncertainty rela- 
tion these momentum space correlations reflect the spa- 
tial and temporal structure of the source from which the 
two particles were emitted. This correspondence forms 
the basis for Bose-Einstein interferometry in nuclear and 



particle physics [Q , a variant of the well-known Hanbury 
Brown-Twiss intensity interferometry in astrophysics || . 

(2) For both pairs of identical and non-identical par- 
ticles additional two-particle correlations can be gener- 
ated by two-particle final state interactions, notably by 
the long-r ang e Coulomb interaction if the particles carry 
charge p-pH|. These final state effects also depend on 
the spatial and temporal distance between the emission 
points of the two particles and thus must contain in- 
formation about the space-time structure of the source. 
While in the last few years we have seen considerable 
progress in our understanding of how to extract from two- 
particle correlations quantitative information on both the 
geometric and the dynamic space-time structure (sizes, 
lifetime and expansion velocities) of the source in the ab- 
sence of final state interactions (FSI) P]l3[], not much is 
known about how to quantitatively correct these methods 
for FSI effects. This is largely due to the lack of a general 
and exact expression which relates the measured correla- 
tion function to the emission function of the source. In 
the absence of FSI (and for chaotic source) this relation 
reads @-£e|§ 



C(q,K) = l± 



If e^ x S(x,K)\ 



f x S(x,K+i)f y S(y,K-l) 



(1) 



where S(x, K) is the single particle Wigner density 
("emission function") of the source, and where 



K 



\{Va + Pb) 
q = Pa - P6 , 



K° = \{E a + E b ), 
q° = E a -E h = a- K/K° , (2) 



V 2 ab - The second term in (|T|) reflects 



with E a ,b = ^Jm* 
the quantum statistical correlations. This expression al- 
lows to expand the correlation function C(q, K) in space- 
time moments of the emission function |l^,[l7| and relate 
the main characteristics of the two-particle correlator, its 
width as a function of the relative momentum q, to the 
second central space-time moment ("r.m.s. widths") of 
the emission function in coordinate space. The factor 
e iq-x m corre i a ti on term in (Q) reflects the assumed 
absence of FSI, i.e. plane wave propagation. 

The goal of the present paper is to generalize rela- 
tion (Q) to the case including two-particle FSI. Existing 
treatments of this problem [0-|^,|l^,[12| differ in their re- 
sults even on the formal level because different types of 
approximations are used already in intermediate stages 
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of the calculation. This, unfortunately, makes a direct 
comparison of these results and a check of the approxi- 
mations essentially impossible. Our aim is to give a rig- 
orous derivation of the generalization of Eq. (|l|), using 
only the following two approximations: 

(1) the source is completely chaotic, i.e. the particles 
are emitted independently, and 

(2) finite multiplicity corrections [6|Jl9[| can be ne- 
glected. 

Both approximations are expected to be good for high 
energy nuclear collisions with large multiplicities, but 
may require further scrutiny in lower multiplicity e + e~ or 
hadron-hadron collisions. On the other hand, high mul- 
tiplicities can result in modifications (screening) of the 
two-particle (Coulomb) potential due to the influence of 
the environment of other charged particles. It was, how- 
ever, shown in [fl8f that such screening effects are small 
and can be neglected because the charge density of the 
environment drops rapidly as a function of the distance 
from the collision center (faster than 1/r 2 ) and the pair 
leaves very quickly the region of high charge density. 



II. THE TWO-PARTICLE CROSS SECTION 

A. General strategy 

The two-particle correlation function C(p a ,Ph) = 
C(q, K) is defined as a ratio between the normalized in- 
variant two-particle coincidence cross section P2(p a ,Pb) 
and the product of the invariant single particle cross sec- 
tions Pi(p a ,&): 



C(q,K) 



-P2(Pa,P&) 
Pl(Pa)Pl(Pb) 



(3) 



where P 2 ( Pa , Pb ) = E a E bl ^-^- b , Pi(p a , 6 ) = E a , b ^. 
Most space will be taken by the calculation of the two- 
particle emission probability P2(p a ,Pb)- Its normaliza- 
tion by the product of single particle probabilities is a 
trivial final step. 

Our strategy is as follows: We start by writing down 
an expression for the two-particle probability amplitude 
Ay(p a ,pf,) for measuring at t — > oo in the detector a 
pair of particles with momenta p Q and p b if the pair has 
been emitted by the source in a two-particle state tp-y. 
The square of this amplitude, averaged properly over the 
distribution of two-particle quantum numbers 7 in the 
source by taking the trace with an appropriate density 
matrix for the source, will give the two-particle proba- 
bility P2(p a ,Pb)- This amplitude, first written down in 
terms of an overlap integral at t — > 00 with an asymp- 
totic two-particle momentum eigenstate in the detector, 
is then rewritten in terms of a corresponding overlap in- 
tegral in the source at the time of freeze-out, by using 
the time evolution operator (including the two-particle 
FSI) to evolve the asymptotic state backwards in time. 



Here a crucial ingredient will be the realization that 
the two-particle FSI can only act while both particles 
are present. If one is emitted earlier than the other, in 
the absence of one-body final state interactions with the 
remainder of the source (e.g. between its charge and 
that of the remaining fireball which we will here neglect) 
it will undergo free time evolution until the second par- 
ticle has also been formed. To implement this requires 
the factorization of the two-particle wave function into a 
product of single-particle states, i.e. the assumption of 
independent emission. This interval of free propagation 
gives rise to an additional phase which contains the time 
structure of the source and which may be important for 
attempts to extract source lifetimes from the two-particle 
correlator. As far as we know, this aspect has not been 
fully discussed in previous FSI studies f^-[T^] . 

Once the two-particle amplitude ^4 7 (p a ,pb) is ex- 
pressed in terms of an overlap integral at the time of par- 
ticle freeze-out, one can calculate the two-particle prob- 
ability by averaging over the quantum numbers 7 and 
the distribution of emission times. Due to the assumed 
factorization of the two-particle wave function the result 
can be expressed in terms of the single particle Wigner 
density of the source at freeze-out S(x, K) which will 
be defined in Eq. ([hj]) . The remainder of the derivation 
consists of an appropriate rewriting of this expression 
which allows to compare it with Eq. (^) and which can 
serve as a starting point for further approximations in 
order to compare with previously published expressions 
in Refs. 

We should stress that our derivation is intrinsically 
non-relativistic. For this to be appropriate we must 
work in a reference frame in which the pair moves non- 
relativistically, i.e. where K«0. Since we are interested 
in the deviation of the correlator C (q, K) from unity for 
small q, the relative motion is then also non-relativistic. 
In this frame the two-particle FSI can be represented by 
an instantaneous potential, and the difference between 
the emission times of the two particles is well-defined. 
The final result can then be evaluated in an arbitrary 
frame by proper Lorentz transformation of all momenta. 
A completely covariant derivation (which should also be 
possible by introducing propagating fields for the FSI) 
has not yet been achieved. 



B. The two-particle momentum amplitude 

Let us consider a two-particle state ip-y emitted by the 
source. Its propagation to the detector is governed by 
the Schrodinger equation 



. i9yj 7 (x a ,x b ,t) 
dt 



ff(x a ,x 6 )^ 7 (x a ,X6,f) , (4) 



where 

P(x a ,x b ) = h(x a ) +h(xi) + V (|x a -x 6 |) , (5) 
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and 



(6) 



The index 7 denotes a complete set of 2-particle quan- 
tum numbers. (In a basis of products of two wave packets 
these could contain the centers X a , Xj, of the wavepack- 
ets of the two particles at their freeze-out times t a , t b , 
respectively.) Eq. (0) is solved by 

^ 7 (x Q , x b , t) = e -<*(*.«)(*-*o) ^ 7 ( Xa , x 6 , t ) (7) 

in terms of the two-particle wave function at some ini- 
tial time to. We will assume that the detector measures 
asymptotic momentum eigenstates, i.e. that it acts by 
projecting the emitted 2-particle state onto 



/out 

h p a ,p b 



' (X a , X 6 , t) = e KP«-Xa-E a t) el ( Pb -x b -E b t) 



— e~ iEt e i (P ci ' Xc '+P f '' X! ') ; (8) 



where E 



E n + Eh 



2K° with E. 



a, 6 



m2 + Pa,fc 



is the total energy of the pair. In this paper we will 
only consider the case of pairs of identical particles, 
m a = rrih = m: correspondingly, the two-particle states 
tp-y in (Q) must be (anti-)symmetrized (see below). The 
measured two-particle momentum amplitude is then 

^ 7 (Pa,Pfc) = 

lim / d 3 x a d 3 x b (j) p u X (x a ,x 6 ,t)?A 7 (x Q ,x b ,i) . (9) 

t— >oo J 

Using the time evolution equation (^) this can be ex- 
pressed in terms of the emitted two-particle wave func- 
tion tpj at earlier times as 

A 1 {p a ,Pb) = lim / d 3 x a d 3 x b ip~(x a ,x b ,to) 

t^oo j 

x exp[-iff(x a ,x 6 )(t o -t)]0° u a * pb (x a ,x () ,t)]* . (10) 

This is correct for all times to > max [t a , t b ] where t atb are 
the freeze-out times for the two particles. Note that by 
partial integration (or inversion of the unitary evolution 
operator) the time evolution operator has been shifted 
from the emitted two-particle state (with arbitrary quan- 
tum numbers 7) to the two-particle momentum eigen- 



statc 



Pa,Pb 



(x a ,x D ,i) which thereby is transformed into 
a distorted wave at time to which includes the effect of 
the FSI in such a way that after evolution from to to 
t = 00 it again becomes a plane wave with momentum 
p a , p 6 and energy E = E a + E b . 



C. The distorted wave 

To further analyze Eq. ( [l0| ) wc introduce center-of- 
mass and relative coordinates according to 



R = -(x Q + x 6 ) , P = 2K = p Q + p b . 
r = x a x fc , q = p a p b , 



such that 



and 



H(x a ,x b )=Ho(R) + H 1 (r), 
£o(*) = -^VL 



(11) 



(12) 



(13) 



with M = 2m and fi — m/2 for identical particles. In 
these coordinates the asymptotic wave function (pf) reads 



^ P6 (x 0) x 6 ,i)=e^ p - R - Bt )e^- 1 



(14) 



In Eq. (|10|) we must evaluate the action of the time evo- 
lution operator on this function. To this end let us sepa- 
rate the total energy E — 2K° into its contributions from 
the cm. and relative motions, E = E cnl + E Te i, where 
E cm = VM 2 + P 2 = 2Vm 2 + K 2 = 2E K « M + P 2 /2A/ 
and E rcl = E-E cia = 2{K°-E K ) » q 2 /(4^) » q 2 /4m 
in a frame where K « 0. The action of Ho(R) on the 
plane wave factor for the cm. motion is then easily eval- 
uated: 



e -iH (R){t -t) e iP R _ g-iVM^+P 2 (t -t) e iP-R _ 

What is left is the time evolution of the wave function 
for the relative motion: 



lim e-^ir^to-t) e ±q-r -i 



E IB it _ -iE„, t 



°0q/2(r). 



(16) 

Here 4>q/2 ( r ) is a solution of the stationary Schrodinger 
equation 



Hi(r) q / 2 (r) = E, eX (j)^/ 2 {r) 



(17) 



(where E Te \ — jfj— ) with asymptotic boundary condi- 
tions 



|r| 



lim q/2 (r) = e^ r 



(18) 



Inserting everything into Eq. ( |l0| ) we find 

A^pi,) = e iEt ° J d 3 rd 3 i? e -' p ' R q/2 (r) 

x V 7 (R+|,R- §,* Q ) . (19) 
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D. Free propagation between emission points 

With Eq. ( |l9| ) we have rewritten the measured two- 
particle momentum space amplitude Ay(p a , p b ) in terms 
of an overlap integral evaluated at the earliest time to 
where both particles are present. Let us, for the mo- 
ment, call this time t a . At this time the first emitted 
particle of the pair has already propagated for a time 
t a — tb if it was emitted at t b < t a . During this time the 
first emitted particle cannot "see" the second particle as 
a separate entity, but only as part of the remaining fire- 
ball. Nevertheless, its charge is there, "hidden" among 
the charges of all the other particles in the fireball. 

This stage of the time evolution is obviously very com- 
plicated. The idealization of letting the two particles in- 
teract via a two-body (Coulomb) potential clearly breaks 
down, since the second particle has not been "formed" 
yet. On the other hand, considering the first emitted 
particle as freely propagating during this time is not nec- 
essarily a good approximation either: there should at 
least be the interaction with charge of the rest of the 
fireball which during this time interval is one unit larger 
than after emission of the second particle. However, con- 
sidering both the interaction between the two particles 
and between each one of them and the remaining fireball 
is a non-trivial quantum mechanical tree-body problem, 
at least until the pair has well separated from the fireball. 

We will therefore study here only the simple approx- 
imation that the first emitted particle propagates freely 
until the second particle freezes out. In order to imple- 
ment this into the formalism we must make the assump- 
tion that the two particles are emitted independently, i.e. 
that at t a (up to symmetrization effects) the two-particle 
wave function -0 7 factorizes: 

ipy(x a ,Xb,t a ) = -^=[ip^ a (x a ,t a )i(} li (xb,ta) 



V2 l 



± V> 7a (x & ,t a )"07 b (x Q ,i Q )] 



(20) 



The indices 7a,7& on the I-particle wave functions now 
label complete sets of f-particle quantum numbers. For 
simplicity of notation we will from now on replace 7 0ib 
by a, b. 

Let us first consider the case that the particle in the 
state ip b w as emitted earlier, at time t b < t a . We can 
then write (see Eq. (||)) 



_ -ih(*.)(t a -t b ) 



(21) 



and hence 
ip-y(x a ,Xb,t a ) - 
0(tq - tb) 
V2 



(22) 



±M^t a )e-^^-^Mx a ,t b ) 



In order to evaluate the action of the free time-evolution 
operator we Fourier decompose ip a , 4>b into momentum 
eigenstates: 



d 3 k a 



■^a(ko,*)e 



i(k a -x- 



k 2 



(23a) 
(23b) 



and similarly for ip b with integration variable k^. One 
then obtains 

0{t a - t b ) f d 3 k a d 3 k b 



^(x a ,x b ,t a ) - ^_ j (2 ^ )3 

x Mk a ,t a ) Mk b , t b ) e -*[».*.+«»(*.-*o+«6t 6 ] 

= i(k a -x a +ki,-Xi,) _j_ e i(k b -x a +k a -Xi,) {24) 

where we have written out in the time factor separately 
the contributions from the Fourier transform and from 
the free time evolution according to (p2|). Note that the 
terms ~ t b cancel between them. 
Inserting ( pi[ ) into ( |l9| ) and defining 



Q = k a + kb 

we arrive at 

4; } (pa,p fc ) 



k= i(k a -k b ) 



0{t a - t b ) f d 3 Q d 3 k 



(25) 



V2 



d 3 rd 3 R 



(2tt) 3 (2tt) 3 
' </ 2 (r)±C q / 2 (r) 

v; a (§+k,i a ) ^(§-k,t b ) . 



x el (Q-P)-R e ik-r e iO(Q,k)t 



(26) 



The superscript (a) in the amplitude should remind us 
that the appearance of the second particle (i.e. the cre- 
ation of the "pair") happened at time t a . In deriving 
(pff ) we used the fact that both <^ q /2(r) and </>_ q /2(r) are 
solutions to ( p^7| ) with the same energy eigenvalue, but 
corresponding to boundary conditions (|I8| ) with opposite 
signs in the exponent. From this it follows immediately 
that 



0-q/a(r) = 0q/2(-r) 



(27) 



With this trick the symmetrization effects on the emit- 
ted two-particle state can be absorbed into the relative 
wavefunction describing the FSI. - In ( p6| ) we also re- 
placed the two-particle quantum number index 7 by the 
pair (ab) of single particle quantum numbers and defined 



(Q,k)=£-w(f + k) -w(§-k 



(28) 



The integration over the cm. coordinate R of the pair 
is trivial and yields (27r) 3 i5 3 (Q — P). After doing the 
integrations over Q and r we get (remember that P = 
2K = p Q + p b ) 



A^(p a ,Pb) = d(t a -tb) 



d 3 k 



e in(K,k)t a l* (a k ) 



(2tt) 3 

x MK + k, t a ) ^b(K - k, t b ) . (29) 
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Here 



E. The two-particle cross-section 



n(K, k) = E a - w(K + k) + £ b - w(K - k) , (30) 
and we defined the symmetrized FSI distorted wave 



1 



*(f.*) = ;^(2> k )±*(-?. k ) 



where 



err e 



ikr 



0q/2(r) 



(31) 



(32) 



is the momentum space representation of the distorted 
relative wave function with asymptotic relative momen- 
tum q/2. 

The opposite situation that the particle in the state 
tp a was emitted earlier at time t a and the "pair" appears 
only later at time tb can be dealt with in a similar way. 
One finds 

A%(p a ,p b ) = e(t b -t a ) J ^ e w) t6$ * (q)k) 

x MK + k, t a ) &(K - k, t b ) . (33) 

The total amplitude is the sum of the contributions from 
the two different time orderings: 



A ab (Pa, Pb) = A$(p a , Pb) + A ab(Pa, P 



(34) 



Let us note for later reference that neglecting the free 
time evolution of the first emitted particle from time t b 
to time t n would have resulted in the ansatz 



V> 7 (x ,X6,t a ) 



0{tq - h) 
V2 



lpa{X-a,t a ) 1pb(Xb,tb) ± %l) a {~X-b,ta)ll)b{Xa,tb) (35) 



instead of ( |22j ) for the time ordering (a), and similarly 
with 6{t a — tb) replaced by 8{tb — t a ) for the time ordering 
(b). This would have resulted in the modified amplitudes 



A^iPa, Pb) = 0(ta - h) e™« /|Wk) 



d 3 k 

x e -i"(K-k)t 6 ^(K - k, t 6 ) , (36a) 

CP.. P>) = - *•) ^ / ^3 ** (f > k ) 

x e-^ (K+k)t " ^ a (K + k,i a ) 

x e -MK-k) t( , ^(K - k, t 6 ) . (36b) 



The two-particle cross section is obtained by averag- 
ing (|34|) and its complex conjugate with the density ma- 
trix defining the source. This density matrix is charac- 
terized by a probability distribution for the two-particle 
quantum numbers (a, b) and by a distribution of emission 
times (t a , tb). We write 

ffeCPa, Pb) = Pab,a'b' A* a , h , (p a , p b ) A ab (p a , Pb) , 

ab,a'b' 

(37) 

and make the ansatz 

Pab,a'b' = Vaa' p{t a ,t a ') V h b> p(t b ,t b >) . (38) 

This ansatz factorizes in such a way that independent 
emission of the two particles is ensured. The summa- 
tion/integration in (B7T) is to be understood as 



^ = ^2 J ' dt a dh dt a < dt b 

%b,a'b' ab,a'b' 



(39) 



According to (]34|) the probability consists of four terms 
which we write as 

Pi (Pa, Pb) = P (aa) + P {bh) + P (ab) + P {ba) ■ (40) 

We shall calculate only the first term p( aa ) explicitly. 
The second term is easily shown to equal the first one 
while the last two terms will be shown to vanish. 
Inserting @ into @ yields 

^(Pa,Pb)=/^|^<l>(§,k')«i>*(|,k) 

x J dt a dt b dt a ,dt v 6(t a - t b ) 6(t a , - t v ) 

x i(n(K,k)t -fi(K,k')t /) 

X V ™' p(ta,t a >)4>a(K + k, t a ) fa (K + k', t a >) 

a a' 

x J2 Ubb ' p{t b ,tb')MK--k,tb)fa(K-k',t v ). (41) 

66' 

Let us introduce new time variables 
1 



Y° = -{h + tv) , y° = t b -t b/: 



(42) 



as well as the inverse Fourier representation (see Eq. (B 



^ a (K + k,X° + ^r fc> (K + k',X -^) = (43a) 



d 3 X d 3 x 

TP a (X + |) e ih(K+k)(X° + |a : )-(K+k).(X+|x)] 

[x - §) e -^( K + k ')( x °-^ n )-( K + k 'Kx-ix)] ( 
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i> b (K - k, y° + £) r b , (K - k', Y° - if) = (43b) 

^ (y + |) e 4^(K-k)(y°+i y )-(K-k).(Y+i y )] 
ifc (Y- |) e -4"( K - k ')(^ -3y°)-( K - k ')-( Y -^)] . 
The phase in ([id]) reduces to 

n(K,k)t„-n(K,k')t B ,= (44) 

A [w(K + k') - u(K + k) + u(K - k') - w(K - k)] 



+x 



2K° 



-£[w(K + k') + w(K + k) + cj(K - k') + cj(K - k)] 
while the ^-functions become 

e (x° - y° + ^) 9 (x° -y°- ^) 

= / dcj 1 -^[x°-y°+|(x°-/)] 



27r w + ie 
/ do/ 1 efa /[ x o_ Y ..i (t o_,o ) ] 
27r w' — ze 



(45) 



We now define the single particle Wigner density 
S(X, K) of the source as 

S(X,K) = J d 4 xe lK x p(x° + ^,X Q -^ 

x v °»> ^ ( x + 1) ^ ( x - f) ■ ( 46 ) 

Using the hermiticity of the density matrix one eas- 
ily shows that S(X, K) real. Then we can combine 
Eqs. ©-© to rewrite © as 



A - fc° - uj(K - k) + it K° - fc ' - w(K - k') - ic 
x / d*A'd*ye i ( fc - fc, ^ jr - y J 



x sfA^ + ^M s(y,if-M^ 



(47a) 



To obtain this expression we shifted the integration vari- 
ables in @ by defining fc° = A - uj - w(K - k), 
fc ' = K° - u/ - co(K - k'). 

The second diagonal term p( bb ) can be calculated sim- 
ilarly: 

p<<*W) = / ? g J |f r *(!.i<') **(!,*) 



A + fc° - uj(K + k) + ic A + fc ' - w(K + k') - it 
x f d 4 Xd 4 Ye^ k '^ x ^ 



x S(X,K+ *±^) s(Y,K-±if 



(47b) 



By relabelling the integration variables, 



X 



Y and k 



fc' 



-/fc'. 



and taking into account that $ (iq, — k) = ±$ (iq, k) 
for bosons and fermions, respectively, one easily checks 
the equality 



- p(bb)l 



(47c) 



The cross terms p( ab ) and p( ba ) represent interference 
between amplitudes of opposite time ordering. They va- 
nish by causality Indeed, each cross term contains under 
the integral the product of two retarded propagators, e.g. 
l/[K°-k°-oj(K-k)+it] and l/[K° + k ' -uj(K+k')-it] 
in p( af> ), which in the complex fc° and k ' planes have 
poles on the same side of the real k° and fc ' axes while 
the corresponding time-energy exponents in the plane 
wave factor have opposite signs. This latter fact means 
that when doing the fc° and fc° integrations, the corre- 
sponding contours must be closed in opposite half planes, 
thus always missing one of the poles. 

Altogether we thus find 



P 2 (p a , Pfc ) = 2P( aa )(p a ,p h ). 



(47d) 



Equations ( J47| ) are the main result of this paper. In the 
following section we discuss it further and study various 
approximations. At this point let us only note that if we 
neglect the phase factor resulting from the free propa- 
gation of the first emitted particle until emission of the 
second one, i.e. use Eq. ( f36| ) for the two-particle momen- 
tum amplitude, we obtain a very similar expression to 
(47a) where only the terms ui(K — k) resp. — lu(K — k') 
in the two energy denominators are missing. The conse- 
quences will be discussed below. 



III. DISCUSSION AND APPROXIMATIONS 



Expression (47a) for the two-particle spectrum is not 
very convenient in practice, because of the poles resulting 
from the two energy denominators. Their contribution 
can, however, be evaluated analytically by doing the k , 
fc ' integrations using residue calculus. For A — Y° > 
one can close the two integration contours in such a way 
that both poles are encircled and contribute. One then 
obtains 



^2(Pa,Pb) = 2 



d 3 k d 3 k' 



$(§,k') S*(f,k) 



(2tt) 3 (2tt) 3 
d 4 X d 4 Y 9(X° - Y°) e *(*-*'M*--r) 

S(X,K+ S (y, A - ^f- 



(48a) 



where 



fc° = K°- lu(K - k) , fc ' = A - lu(K - k') (48b) 



G 



if the free time evolution between emission points is in- 
cluded, while 



k° = k ' = K° , 



3c) 



if it is neglected. 



A. Wigner representation 

Equation (|47|) can be written as a folding relation 
between Wigner densities. We define p = \{k + k'), 
Q = k — k' , and 



Xq/2 (p)^$(§,p) 



K° + lu{p-K) 



(49a) 



for the case where the propagation between emission 
points is included, or 



X q/2 (p) = *(§, p) -p- 



1 



K° 



(49b) 



if free propagation between emission points is neglected. 
The two-particle cross section ( |47j ) then reads 



f 2 (Pa,Pf> 



with 



d 4 p 



d 4 X d 4 Y S(X, K + p) 



(2tt)4 

x W n/2 (X-Y,p)S{Y,K-p), 



(50) 



di Q „-* 



Q-(X-Y) 



xx q/2 (p+f) X q/2 (p-§) ■ (51) 

This function W can be interpreted as the Wigner density 
associated with the distorted wave describing the final 
state interactions. It describes the power of the FSI to 
"push" two particles, which were originally emitted with 
momenta p' a — K + p and p b = K — p, to the observed 
values p a and p (resp. K and q). It is determined by the 
"modified distorted waves" \- Please note that the rep- 
resentation (pOh is generic; it was previously derived by 
Pratt in Eq. (2.8) of Ref. with quite different meth- 
ods. Different approximations in dealing with the propa- 
gation of the first emitted particle between the emission 
points only result in different "modification factors" as- 
sociated wi th th e distorted waves <&(t|,p). Equations 
(49a) and ( 49b ) are two specific examples: we expect 
that different assumptions about what happens to the 
first particle while it waits for the appearance of the sec- 
ond one can be similarly included in Eq. ( |5l| ) by simply 
changing the "modification factor" in the definition of 

Xq/2(P)< 



B. Free particle limit 

To discuss the limiting case of no fi nal state interac- 
tions it is best to start from Eq. (|48a| ). In this case the 
wave functions </>±q/2( r ) in jl7| ) describing the relative 
motion of the two particles on their way to the detector 
become plane waves, with the momentum space repre- 
sentation 



5k| = 



V2 



5^ (f 



k)±<5( 3 )(§ + k) 



(52) 



The two "^-functions in (48a) result in four terms two of 
which correspond to ^(k+k 7 ) = ±q, k— k' = 0, while the 
other two terms have |(k + k') = 0, k — k' = ±q. The 
corresponding energy values depend on whether or not 
we include free propagation between the emission points. 
If we include it, the energies are + k ') = ±q°, 

k° - k ' = for the first two terms and |(fc° + fe ') = 0, 
k° - k ' = ±q° for the last two terms (q° = E a - E b ). 
The two-particle cross section becomes 



^ 2 (Pa,P 6 )= J d*Xd 4 Y6(X -Y°) 



S(X,K+1) S(Y,K-$) 
S(X,K-i) S(Y,K+i) 



± 



,iq-(X-Y) 



,-iq-(X-Y) 



S(X, K)S{Y 7 K) 



(53) 



After relabelling X Y in the second and fourth term 
and using 0(X°-Y°) + 6(Y°-X°) = 1 we obtain 



^ 2 (Pa,Pfc)= / d 4 Xd 4 Y [s(X,K+i) S(Y,K- 



± e iq<x-Y) s(X,K)S(Y,K) 



(54) 



After normalization with the single particle spectra we 
thus recover the correct expression ([j ) for the correlator 
in the limit of vanishing final state interactions. 

If the free propagation between emission points is ne- 



glected, according to (48c) all four terms have the same 
energies ^(k° + k° ) — K°, k° — k° = 0, and it is obvi- 
ous that it is not possible to recover an expression with 
any close similarity to (p)4|). The phases resulting from 
the time evolution of the first particle until the appear- 
ance of the second one are thus crucial for the correct 
free-particle limit. One may be able to modify the prop- 
agation between emission points, but not to drop it com- 
pletely. The correct free particle limit requires the ap- 
pearance of a phase ~ (t a — tj,). 



C. The smoothness approximation 

For practical applications one would like to know how 
to combine this formalism with classical event genera- 
tors which determine the emission functions S(x, K) as 
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a set of last interaction points in the kinetic evolution 
the many-particle collision system. For particles with- 
out FSI it was shown in j| that the symmetrized two 
particle cross section (|^) can be computed by sam- 
pling the generated particles at the on-shcll momenta 
(p a ,E a ), (pb,Eb) for the direct term, and at the on- 
shell value (K, Ek) for the exchange term (with both 
selected particles having the same on-shell momentum 
(K, Ek))). The pairs selected for the exchange term 
are then multiplied with a "Bose-Einstein weight factor" 
cos q ■ (xi — Xj) where q is the relative momentum value 
at which we wish to know the correlator (not that of the 
selected pairs which have pi = pj = K), and 
the space-time distance between the particles in the se- 
lected pairs. This doe s no longer work in the presence 
of FSI. Equations (48a), ( |5C| ) show that now we need to 
sample the emission function S(x,p) at all possible and 
in general off-shell values of the momentum p. The parti- 
cles become only on-shell when they reach the detector, 
by virtue of t he F SI. It is thus basically impossible to 
simulate Eqs. (48a), ( (50| ) with classical event generators 
unless one imposes a further approximation which essen- 
tially eliminates all off-shell effects. 

This step is known as the smoothness approxima- 
tion in which one assumes that S(x, K ± p) has a suf- 
ficiently weak momentum dependence that, over the p- 
range where W^^iX — Y,p) in Eq. ( |50| ) is non- vanishing, 
S(x, K±p) can be replaced by either S(x, K) or S(x, K± 
p), depending on what seems more convenient. The do- 
main of support of the function W q /2 (X — Y, p) is a mea- 
sure of the ability of the FSI to change the particle mo- 
menta after freeze-out. For our purposes the smoothness 
approximation can be considered reliable as long as the 
typical shift in momentum in the FSI is small compared 
to the q-range over which the correlator shows interest- 
ing structure. Note that for massive particles the FSI 
cause mostly a change of the spatial momentum while 
the corresponding energy transfer is very small. This is 
the basic reason why off-shell effects are small and why 
the smoothness approximation works. 

Our implementation of the smoothness approxima- 
tion differs in a crucial detail from previous approaches 
B'illllJUl- By formulating the approximation concisely 
and implementing it systematically we obtain an expres- 
sion which shows a strong similarity to the free particle 
case in that it evaluates for the direct and exchange terms 
the source function at different momenta. In this way we 
automatically avoid some of the possible pathologies in 
the behaviour of the cor relator which may arise 

from the more traditional versions of the smoothness ap- 
proximation. 



x g-^KK-kK^K-k')]^ e -i(k-k')-y 



We start from Eq. (48) which we rewrite as 



P2(p„,Pl) = 2 fd'xd'yeiy' 1 ) 



x s(x + Z,K+*&-) s(x-hK- . (55) 



Let us expand the frequency difference in the temporal 
phase factor for small values of k, k', keeping terms up 
to second order: 



[v(K - k) - iv{K - k')] 



(56) 



Using (|3l|), the two-particle spectrum thus becomes 
ft(p«,Pi) = / d i xd 4 y6{y <) ) 



d 3 k d?k' .(k-kQ.fy-^K-^y) 

(2tt) 3 (2tt) 3 

(f,k) 0(f,k') ±^(f,k) 0(-f,k') 



±< 



-f,k)0(§,k') 



(-f.k) 



x s{x + l,K+^f\ s(x-%,K-^f 



(57) 



The functions </> describe the probability amplitude of 
finding in the FSI distorted wave with asymptotic relative 
momentum ±q/2 a plane wave with momentum k resp. 
k'. These functions are peaked around k, k' = ±q/2, 
and the peaking is the stronger the weaker the final state 
interactions are. The four terms in the square bracket in 
@ thus peak at (k+k')/2 = q/2, 0, 0, and -q/2, respec- 
tively. (One easily checks that these equations hold for 
the corresponding 4-vectors.) We will use this property 
by replacing (k + k')/2 in the phase factor and in the ar- 
guments of the emission functions by the corresponding 
peak locations. This allows to pull the emission func- 
tions outside the k, k' integrations; it also removes the 
quadratic dependence of the phase factor on k, k' in fa- 
vor of a linear one, turning the integrations over k and 
k' into normal Fourier integrals. The Fourier integrals 
can be performed, giving the corresponding relative wave 
functions in coordinate space: 



P 2 (Pa,p 6 ) = J d 4 y6(y°) 



(58) 



q /2(y-v 6 y°)| 2 / d 4 xS (x + %,p a ) S(x-%,p b ) 



±2 Re (^ /2 (y-vy°)0_ q/2 (y-vj; ) 

x J d 4 xS(x + %,K) S(x-%,K) 

+ |0-q/2(y-V a 2/°)| 2 / d A xS(x+ |,p 6 ) S(x-%,p a ) 



Here we defined the three velocities 



K 



Pa I"/ . r , 

= , v a = — , v 6 = — (59) 



Pb 



E K 



E k 



E k 
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associated with the observed particle momenta p a , p , 
and their average K. 

Using the relations ( |27| ) one shows that in the middle 
term in ( fisl) the factor containing the FSI relative wave 
functions is even under a simultaneous sign change of y 
and y° . The product of emission functions in this term is 
even under the same sign change, too, which means that 
in this term we can replace 0(y°) by ^[8(y°) + 9(— y )] = 
i. The first and last term in ( |58| ) do not have such a 
symmetry, but can be combined by substituting y — > —y 
in the last term to yield the final expression 

P2(Pa,Ph) = [ dPxdPyS (x+ \,Pa) S (x - 

x %°) |0 q/2 (y-v 6 y o )| 2 + d(-y°) K/ 2 (y-v a2/ °)| 



± J d i xd 4 yS{x + ^,K) S(x-%,K) 

x C q /2(y- v 2/ o )0q/2(y-vt/°). 



(60) 



Once again one checks that in the limit of vanishing FSI 
(i.e. by replacing the functions <fi by plane waves), and 
using the approximation K° w E% in accordance with 
the approximation made in Eq. (p6|), the correct expres- 
sion (|54|) is recovered. 

It is worthwhile to discuss the physical meaning of the 
arguments at which the FSI distorted waves <j> in J6C| ) 
must be evaluated. The sketch presented in Fig. [j] shows 
the relevant possibilities for the first term in (|6(]) (the "di- 
rect term"). In this term the two emission functions are 
evaluated at the observed momenta p a ,Pb, corresponding 
to particle velocities v a , Vb. Fig. [j] shows that in each case 
(i.e. for both possible time orderings between the two 
emission points) the argument of the FSI distorted wave 
cj> corresponds to the spatial distance of the two particles 
at the time when the second particle freezes out, i.e. when 
the pair first exists. The different velocities which arise 
for the two different time orderings reflect the velocity 
of the earlier emitted particle in each case. In the sec- 
ond term of ( |60| ) (the "exchange term" ) the two emission 
functions are evaluated at the same momentum, namely 
the average pair momentum K, and correspondingly for 
both time orderings the earlier emitted particle has the 
velocity v corresponding to this momentum. Again the 
argument of the FSI distorted wave is the spatial distance 
between the two particles at the time of emission of the 
second particle. 



D. Implementation in event generators 

Equation (^) can be easily implemented into event 
generators, following essentially the same procedure as 
given in Ref. j|: 

For the direct term one selects all pairs with 
Pi = Pa , Pj = Pb within a given numerical accuracy (bin 
width) which is essentially dictated by event statistics. 



Each pair is multiplied with a weight given by the corre- 
sponding probability density |0 q /2| 2 of the FSI distorted 
wave. The latter must be evaluated in a frame in which 
the pair moves non-relativistically, best in the pair rest 
frame where K = (p a + p&)/2 = 0. (Then Ek = m, 
and the velocities ( |59| ) reduce to their usual nonrelativis- 
tic definition.) From the space-time coordinates X4, Xj of 
the pair in the event generator frame one calculates the 
distance y*j between the two particles in the pair rest 
frame at the time when the second particle is produced. 
One then computes |</> q . /2(y*j) | 2 and weights the selected 
pair with this number. In this expression q* is the 
spatial relative momentum between the two particles in 
the pair rest frame which must be computed from p a , pt 
in the event generator frame. (We remark that in the 
absence of FSI, the corresponding weight is simply 1.) 
The complete direct term is obtained by summing over 
all such pairs. 

For the exchange term, the selection of pairs and 
weights is a little less intuitive ||: One selects all pairs 
(i,j) with pi — pj — K (i.e. on-shell particles (!) with 
Pi = p j = K and £7, = Ej = Ek), again within the same 
numerical accuracy (bin width) as above. From the pro- 
duction coordinates X4 , Xj one again computes the spatial 
distance y*^ between the two particles in the pair at the 
time of emission of the second one, in the pair rest frame 
K = 0. This distance is used to compute the weight 

't'-q' /2(.ytj) < l > q* /2(ytj) f° r this P au "- The v& l ue 01 h ere 
is the same as above in the direct term, i.e. it is computed 
from p a and p by transforming into the pair rest frame, 
not from pi = pj = K. (Without FSI, the corresponding 
weight would be cos(q* - y*j) ||-) The full exchange term 
is obtained by summing over all such pairs. 

Finally one must normalize the correlator by the prod- 
uct of single particle spectra, 



Pi(p a ) Pi(Pb) = J d 4 xS(x,p a ) J d 4 yS(y, Pb ). 



(61) 



This normalization is best obtained from the pairs se- 
lected for the direct term above by multiplying them with 
unit weights. 

One may object to the use of event generators for the 
emission function because they fix particle momenta and 
coordinates simultaneously and thus violate the uncer- 
tainty principle. It was shown in || how to generate 
from an event generator a quantum mechanically con- 
sistent Wigner density S(x,p) by folding the event gen- 
erator output with minimum uncertainty wave packets. 
The corresponding quantum mechanically consistent al- 
gorithm for computing single- and two-particle spectra 
given in || is easily generalized to include FSI effects, 
by simply replacing the factors 1 and cos(q* • y*.- ) in the 
direct and exchange terms, respectively, by the correct 
FSI weights as discussed above. 







E. Comparison to previously published expressions 



F. Nonidentical particles 



In previous treatments gimjH] the smoothness ap- 
proximation was implemented in a different way: the 
smoothness of S(x,p) as a function of momentum was 
used to replace the momentum arguments of the two 
emission functions in both the direct and the exchange 
term of Eqs. ([54]), (60) in exactly the same way, namely 
either by p a in the first and by pb in the second, or by K 
in both emission functions. The first of these two alterna- 
tives can, for sources with very strong x-p-correlations, 
lead to correlators with pathological behaviour, as dis- 
cussed in . The second alternative was exploited 
in the FSI studies presented in Refs. . In this case 
both the direct and the exchange term in (|60|) involve the 
same combination of emission functions, namely 



D(y,K) = J d A xS{x + ^,K) S(x-%,K) 



(62) 



D(y,K) is the "relative distance distribution" of the 
source, i.e. the distribution of relative space-time dis- 
tances y between the particles in emitted pairs with 
momentum K. Obviously D is an even function of y, 
D(y, K) = D(—y,K). Also, since in the direct term of 
( |60| ) the emission function arguments p a ,Pb are now both 
replaced by K, the velocities v a ,Vb in the corresponding 
FSI distorted wave factors must also both be replaced 
by v. The two time orderings can then be combined by 



using 



9(y°) + 9(-y 



1, and we obtain 



P 2 (p a ,p 6 ) =P 2 (q,K) 



d 4 yw(y;q,K)D(y,K), 



with the "weight function" 
w(y;q,K) = |$ q / 2 (y-V2/ )| 



(63) 



(64) 



|<£ q/2 (y-vy o )±0 q/2 (-(y-vy o ))| 



[Note that if P 2 is approximated as in (^3|) one should 
also for consistency use the corresponding approximation 
Pi(p a )Pi(pb) = [Pi (KYI for the normalization of the 
correlator.] Eqs. (j63|), ( |64| ) are identical with Eq. (2.11) 
in Ref. |11|. In the pair rest system (v = = K) P 2 



reduces to 

P2(Pa,Ph) 



For correlations introduced by FSI between non-iden- 
tical pairs one has to take into account two differences: 

(1) the two masses m a and m\, are unequal, and 

(2) the two-particle wave-functions are not symme- 
trized resp. antisymmetrized. 

The first modification requires the introduction of corre- 
spondingly modified center-of-mass and relative coordi- 
nates in several steps of the derivation given in Sec. 0. 
Still, the problem remains non-relativistic in any system 
in which the pair is at rest or moves non-relativistically, 
and the rest of the derivation goes through as before. The 
final result is again given by Eq. ( poj ) ; the only difference 
is that now <£> q / 2 (y) is just the distorted wave for the 
relative motion including FSI, not its symmetrized form 
( pl| ) . As a result of the mass asymmetry, the weight func- 
tion w(y; q, K) even in the pair rest system will no longer 
be symmetric in y. Its asymmetry will be reflected in an 
asymmetry of the correlator under sign change of q which 
can be used to extract information about the average dis- 
tance between the effective sources from which the two 
particle species originate ]24|]. 



G. Coulomb final state interactions 

For completeness and to correct a few confusing ty- 
pographical errors in Ref. @| we discuss explicitly the 
distorted wave for Coulomb final state interacions. In 



this case the solution of the Schrodinger Eq. (17) is 

0±q/a(r) = T(l + irj)e-^ e ± ^' r P(-^ ; , (66) 

where (q = |q|, r = |r|) 

afi am 



V = 



q/2 q 
is the Sommerfeld parameter, and 

Z T = | (l r T q • r ) = \qt(1 T cosfl) . 



(67) 



(68) 



being the angle between q and r. F(—im 1; iz) is a con- 
fluent hypergeometric function with the series expansion 



d 3 y |$ q/2 (y)| dy°D(y,K). (65) F{-i m l;iz) = ^ 



T{n-irj) {iz) n 

J r(-«7) W 



This expression was first written down by Koonin ]22| ] 
and has recently been used in [P3| . Note that on the 
r.h.s. of ( p5| ) everything must be evaluated in the pair 
rest system with K = 0. In this version of the smooth- 
ness approximation the two-particle spectrum P 2 (q, K) 
is thus given by the time integrated relative distance dis- 
tribution in the pair rest system defined by K, weighted 
with the probability density of the outgoing distorted 
wave with relative momentum q at the spatial relative 
distance at which the pair was created. 



1 — ii] — if?)(2 — if]) ■ • • (n — 1 — irf) 



{iz) n 



(69) 



For small 77 1 this can be approximated by 

n=l 

= 1 -ir)(Ci(z) +iSi(z) -Inz -7) . (70) 
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This is good for pions with q ^> 1 MeV and for protons 
with g> 7 MeV. The square of the normalization factor 
in (66) yields the well-known Gamov factor 



G(q) = 



2nr] 



>2ttJ7 _ ^ 



r(i + «?)| 2 



(71) 



With these expressions the two-particle cross section ( pOj ) 
takes the form 



P 2 (q,K) = G(g) |l 1 (q,K)±J 2 (q,K) 
where the direct term is given by 
/i(q,K)= / d 4 y 



(72a) 



%°) \F(- im l; iZ b _(y ))\ z 

+e(-y°) \F (- i77 ;l;zz°(y ))| 2 ] (72b) 
x f d 4 xS(x + %,K+%) S(x-*,K-Z) 
while the exchange term is 

/ 2 (q,K) = J d 4 ye l «<y- v y°^D(y,K) (72c) 
x F(ir,;l;-iz+(y )) F {-i m l:iz_(y )) . 



Here D(y,K) is the relative distance distribution (62), 
and z_ (y°), z"(y°), z±(y°) are determined with the 
help of ( |68| ) by substituting r — > y— v^y , y— v a y°, and 
y— vy°, respectively. Note that the integrals 7x,2 should 
be evaluated in_(or close to) the frame where K = v = 0. 
Equations (||) correct Eq. (3.11) in Ref. @]. 



IV. CONCLUSIONS 

Let us shortly summarize our main results: 
We studied the effect of two-body final state interac- 
tions on the two-particle coincidence cross section, both 
for pairs of identical and of non-identical particles. We 
used only two approximations: (1) the source is chaotic, 
and (2) the total particle multiplicity is large. The first of 
these two assumptions is crucial since it allows to factor- 
ize the two-particle Wigner density of the source and ex- 
press the two-particle cross section in terms of the single- 
particle Wigner density, i.e. the emission function S(x,p) 
of the source. 

From these two assumptions we derived, without fur- 
ther approximation, the general expression (|4^) for the 
two-particle cross section in terms of the emission func- 
tion S(x,p) and the FSI distorted relative wave func- 
tions in momentum space, $ (?f,k). This expression can 
be rewritten in the generic Wigner representation (50). 
In spite of its generality, this expression is of limited 
practical usefulness because it involves the emission func- 
tion at arbitrary of-shell momenta which is usually not 
known. Since for massive particles the FSI shifts, how- 
ever, mostly the spatial momenta of the particles while 



the corresponding energy transfer is very small, it is in 
practice possible to replace the off-shell momenta in the 
emission function by on-shell values. We do this sys- 
tematically in the framework of the so-cal led sm oothness 
approximation which we discussed in Sec. [II C. Here our 
treatment differs from previously published ones, and our 
formul a (pG ) thus improves upon known results. 

Eq. (|6fJ) involves only on-shell momenta and can thus 
be implemented i n clas sical event generators. This was 
discussed in Sec. [IIP. In contrast to previously pub- 



lished expressions, our result exhibits the same asymme- 
try in the momentum arguments of the emission func- 
tions between the direct and exchange terms as expres- 
sion d54| ) for the free case to which it correctly reduces 
when the FSI are switched off. 

This asymmetry between the direct and exchange 
terms is eliminated if one implements the more stringent 
version of the smoothness approximation used by Pratt 
and collaborators which replaces the momentum argu- 
ment of the emission functio n by t he average pair momen- 
tum K everywhere. In Sec. Ill E we showed how in this 
way the well-known Koonin-Pratt formula, Eqs. (|63|)- 
(pq), is recovered. 

The differences between the "correct" asymmetric form 
and the "approximate" symmetric Koonin-Pratt form 
were recently extensively investigated for the case of 
free particle propagation without FSI after freeze-out. 
They were found to be potentially severe for sources with 
strong ir-p-correlations (ejj. sources which feature strong 
collective expansion) H,^0|, but the problem appears to 
be less serious if the source is sufficiently large [^lj . While 
the existence of the "correct" asymmetric expression ( |60| ) 
(for Coulomb FSI an explicit expression was given in 
(|72])) now eliminates the need for using the approximate 
and somewhat uncertain Koonin-Pratt formula, it would 
still be nice to have a quantitative feeling for the error 
margin associated with the Koonin-Pratt approximation. 
A systematic numerical study of this question which in- 
cludes FSI effects is under way. 
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